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[57] ABSTRACT 

A computer implemented tool forecasts the spot price of 
electric power in a deregulated market and the amounts of 
power that may be traded in this market. Using generating 
capacities of multiple utilities, price functions, weather 
forecasts, and transmission variables, the computer imple- 
mented tool makes these forecasts at different delivery 
points, providing the decision maker with probabilistic 
distributions for spot prices and trading. 
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POWER TRADING AND FORECASTING unbundling the electric power into generation and transmis- 

TOOL sion. In the future, one will pay a production cost and a 

transmission fee. There will be several power suppliers from 

CROSS REFERENCE TO RELATED whom electric power may be purchased. Suppliers may have 

APPLICATION 5 different pricing mechanisms. For instance, there might be a 

discount for using power off-peak periods or for signing a 

This application is related in subject matter to U.S. patent long-term contract with the supplier. Power producers will 

application Ser. No. 08/869,561 (IBM Docket YO997-032) compete with each other to minimize their costs so that they 

filed Jun. 5, 1997, by Samer Takriti for "RISK MANAGE- can sell their product to more customers and maximize their 

MENT SYSTEM FOR ELECTRIC UTILITIES" and profit. 

assigned to a common assignee. The disclosure of applica- The hope is that deregulation will result in cheaper prices 

tion Ser. No. 08/869,561 is incorporated herein by reference. and play a part in improving the economy by encouraging 

investments in electric utilities. The size of the electric 

DESCRIFnON industry is expected to grow after deregulation as was the 

. ^rs^^^r xvrr^ ^ ™ TT . ¥VT , ^vt^t^v, c case with the telecommunications industry. The telecommu- 

BACKGROUND OF THE INVENTION « nications Kyemt shot up ^ m m bmion tQ 

1. Field of the Invention $170 billion within ten years of deregulation. 

The present invention generally relates to computer Apitfall of deregulation is that the load on a utility system 

implemented forecasting tools and, more particularly, to a ™ becoming increasingly unpredictable. The reason is that 

tool for forecasting the spot-market prices of electric power 20 tradin g transactions can change the load pattern signifi- 

and the trading transactions at different delivery points. cantlv ; For example, some utilities may sell more than 30% 

/■> t% i j r» • i* of their power generation to other utilities on certain days. 

2. Background Description „ , F , & . . 4 . , 4 - . f , /, 

° Demand and supply in the market are functions of volatile 

Electricity is an essential part of our lives. If we took a electricity prices which in turn depend on highly unpredict- 

moment to think of the services that would not exist if able elements such as reg i 0 nal weather conditions and fuel 

electricity were not available, we would be amazed. Almost prices 

every aspect of our modern lives involves electric power, t*k^« „™ „ n o 

r . , „ .... . . , f Inere are two components tor any competitive environ- 

trom light bulbs and television sets to hospitals and auto- „„♦ «u , • i a 1™ • i t *u pi * • *u 

. . , . A ,, , , , . ment: physical and commercial. In the case of electricity, the 

mobile industries. Although we are used to having power ^u.,„- * i ♦ „ • . * *u e c . • 

* . i j . & / , . physical component already exists in the form of transmis- 

whenever we need it, the processes and systems involved in sioQ Unes and aU consumers have h ysical access to the 

dehvenng electricity require careful planning and sophisti- market ^ coamtchd component mnskls of the ability of 

cated mathematical and forecasting models. buyers aQd suppHers tQ execute ^ ovef ^ phy J cal 

Lately, due to rising costs of energy and discrepancies in network. One should also ensure that anyone who is inter- 
its price among different regions of the country, the legal e sted m joining this market has non-discriminatory access to 
framework surrounding the electric-power industry has 35 the physical network. In this environment, suppliers and 
changed. This change has opened the door for utilities to consumers have the freedom to contract with each other 
compete with each other and against independent suppliers unc j er negotiated terms and conditions. To ensure the deliv- 
regardless of their geographic location. Although this eryof the power traded, they need a mechanism to ensure the 
change will benefit the consumer, utilities are going to face physical delivery of the electricity. Also, they may need to 
a highly unpredictable market and will need to make tough 4Q rely on a sp 0t market to buy or ^ additional electricity, 
decisions regarding power generation and delivery. Although the new environment encourages competition, 

The power industry is going through deregulation. The power transmission will remain regulated for the time being 
current picture of a single utility controlling the market in a m or d er to maintain a reliable system. The transmission lines 
specific region will soon disappear. Instead, there will be m ea ch state or region will be controlled by an independent 
power producers who sell their product to a power pool; and 45 entity called Independent System Operator or ISO. One of 
power suppliers who will buy power from this pool and in the ISO responsibilities is to settle financially with the 
turn sell it to their customers. Although the full picture of the part i es involved in transmitting electric power. The trans- 
power industry after deregulation is not yet known, it is clear mission cost depends on the proximity of the supplier and 
that utilities need to prepare themselves for an open market the congestion of the transmission lines as well as other 
in which buying and selling power are to be considered 50 operational factors. To maintain a reliable system, ISO 
when scheduling their generating units. announces, twenty-four hours in advance, the load forecast 

The main reason behind deregulation is to reduce the high on the system and asks interested suppliers to submit bids; 

price of electric energy. Initial steps towards deregulation i.e., amount of power to be delivered at each point and 

were taken in 1978 with the passage of the Public Utilities pricing information. The load forecast is given for each 

Regulatory Policy Act. This act encouraged nonutility gen- 55 delivery point within the physical system, 

eration and required utilities to buy power from independent Bids are submitted to the Power Exchange or PX which 

generators. The Energy Policy Act of 1992 took deregulation then holds an auction to determine which suppliers to buy 

a step further by mandating open access to the transmission power from. For each delivery point, the PX chooses the 

system for wholesalers. least expensive source of power, then the second least 

Currently, electricity is sold as a service that is delivered 60 expensive, and so on. The process is repeated until the 

to specified points. For example, each one of us expects to forecasted load on the system is met. Note that the previous 

receive electric power via a meter outside the house. We pay process ensures that the customers' need for electricity is 

for this service regardless of its producer or which power satisfied at a minimal cost. The bidding cost of the last 

lines it followed. That is, an electricity bill indicates the total supplier chosen in the bidding process is the spot-market 

usage of electricity in kilowatt hours (KWH) and the service 65 price for power at this delivery point. The PX performs this 

price per KWH without incorporating any other details into process for each delivery point in the system and at each 

the pricing scheme. Deregulation is changing this picture by hour of the day. 
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After agreeing on the production schedule, the ISO stud- different amounts of electric power. The likelihood of a 

ies the schedule to make sure that it is feasible; i.e., it does transaction is given in the last column, 

not exceed the capacity of the transmission lines. Transmis- Clearly, one would like to use all available historical, 

sion charges are a function of the congestion of the power competitive, and weather information to produce the most 

lines. These charges are calculated by ISO and billed to the 5 accurate output. However, the data available for different 

different parties involved in trades. utilities may differ from one region to another. For example, 

Note that suppliers who submit a bid with a high price so™ utilities are required, by the power pool, to post the 

may end up not selling any of their production. On the other status of their generating units by the end of the week. These 

hand, selling power at a low price may not create enough requirements differ significantly between power pools. Also, 

revenue for a generator. Given that no one knows in advance 10 data that are currently public, such as hourly loads, may not 

the amount of power that competitors may bid for, the De available in the near future, 
electric-power market will become more uncertain and 
risky. The hope is that competition will force utilities to 

minimize their production costs which will, in turn, benefit The foregoing and other objects, aspects and advantages 

consumers and the economy in general. The previous model *5 win be better understood from the following detailed 

also encourages individuals to invest money in the electric- description of a preferred embodiment of the invention with 

power market. For instance, if the charges across a certain reference to the drawings, in which: 

power line -are relatively high due to congestion, an investor FIGS 1A and 1B are block diagrams uiustrating the 

could build an alternative power line in the hope of ben- relation of the power trading and forecasting tool according 

efitting from the high-traffic in that region. The result, of 20 io ±& mvemion tQ {{s m and 

course, is a reduction in the transmission cost in that region _ T ^, . . ,. .„ - 

and an increase in the reliability of the system. FIG -. 2 * » flow . ^ »U«f ^8 the process of 

approximating the price function of a utility by a quadratic 

SUMMARY OF THE INVENTION function; 

It is therefore an object of the invention to provide a tool 25 , FI ,? S - 3A and 3B , are demonstrating how the cost 

for forecasting the spot price of electric power in a deregu- ^ n ^ 0nof 3 generat0r Can be a PP r0Xlnlated US1D S a llDear 

lated market and the amounts of power that may be traded nc l0n ' 

in this market FIG. 4 is a graph showing an approximate cost function 

The computer implemented process according to the 3Q o£ a utmty constnicted by usmg tte linear approximations of 

invention forecasts the spot-market prices and the trading * ' 

transactions at different delivery points, providing the deci- FIG - 5 is a flow diagram of the simulation used in the 

sion maker with probabilistic distributions for spot prices process according to the invention; 

and trading so that he or she can manage risk efficiently. The FIG. 6 is a flow diagram of a process that generates a 

user of the process can be (1) a utility or an independent 35 random sample; and 

power supplier who is interested in using the trading and FIG. 7 is a flow diagram showing the bucketing process 

spot market information to decide which generating units to used to create a probabilistic distribution using the output of 

commit, how to price his or her power, and the level of simulation runs, 
generation needed; or (2) a power broker who is interested 

in buying and selling power (but not producing) in order to ^ DETAILED DESCRIPTION OF A PREFERRED 

maximize his/her revenue. The computer implemented pro- EMBODIMENT OF THE INVENTION 

cess uses the weather forecast in addition to other market Referring now to the drawings, and more particularly to 

information to provide probabilistic distributions for the FIGS. 1 A and IB, there are shown in block diagram form the 

power traded and spot-market prices. What distinguishes our relat ionship of the power trading and forecasting tool 

tool from others is that it provides a distribution instead of 45 according t0 the invention and its inputs and outputs. The 

a smgle estimate for the spot price and the power traded. For lool can nave differem sels of input parameterSj and two are 

example, the computer implemented process may provide desc rib e d here by way of example. The first example is 

the user with Table 1 so that he or she can hedge against shown in FIG 1A {n this ^ we are given tne historical 

uncertainty in the market. data of prices and amounts of po Wer sold by each 

so utility. The tool 11 receives inputs from a plurality of utilities 

TABLE 1 i2 a to 12 n . These inputs include for each utility prices at 

which it sold power in the past and the corresponding 
amounts of power sold. In addition, the tool 11 receives as 
an input the transmission costs and the maximum transmis- 
55 sion capacity for different fines in the transmission network 
13. The output of tool 11 is the probabilistic distribution for 
trades and spot prices 14. 

The second set of input parameters is illustrated in FIG. 
IB. They include the operating cost curves of the different 
60 generators 12 i to 12„, estimates of the marginal profits, and 
electric load forecasts. The output is similar to that of FIG. 
1A. The inputs and outputs of the tool 11 are described in 
more detail below. 
System's Inputs 

The first column represents the spot market price. The 65 As we mentioned earlier, the final picture of the deregu- 
second column contains the amount of power that may be lated market is not well defined. Hence, the input 
traded. Note that for the same spot price, one may buy or sell parameters, data available about utilities and their 



$ 


MWH 


Pr, 


12.00 


-50 


10 


12.00 


-40 


5 


12.00 


-30 


5 


12.50 


-40 


15 


12.50 


-30 


10 


13.00 


-30 


5 


13.00 


-20 


10 


17.50 


40 


5 


17.50 


50 


5 


18.00 


50 


15 


18.00 


60 


10 


18.00 


70 


5 



06/22/2004, EAST Version: 1.4.1 



''' ■ .2 0) 



5,974,403 

5 6 

generators, are not exactly known at this time. To avoid this operating environment, another might be the cost observed 

difficulty, we describe the input parameters in general terms when natural gas prices are high, a third may be the price 

so that our tool can handle any electric power market. The function when there is a heat wave in the area, and so on. 

input to our system consists of (1) the generating capacity One can use as many scenarios as needed so that uncertainty 

and price functions, f,- /, of utilities competing in the market, 5 in f is approximated reasonably. 

(2) the load forecast 'in the regions of interest, and (3) the Given a set of observations, \ j-1, . . . , J ( /, of the 

transmission costs, c (iJ) /, charged by the ISO to transmit P rices of utihtv 1 al time 1 under certain conditions s, we find 

electricity from node i to node j at time period t and the lhe P ncin S P rofile that bcst * fits this utUit y usin & a least 

transmission capacity, u (iJ) Note that all of our notations *P m a PP roach - Herc 15 a pseudocode of this process, 

contain the superscript s. This is the scenario index, a 10 Initialization. Get a set of paces, XJf, offered by utility 

concept that is described in some detail in the above refer- \ at time , 1 under s u ^ x «v * £ e ™ ou ? of 

enced application Ser. No. 08/869,561 and which shall be SSs ° 

clarified through the following discussion. Here is a detailed _ P nce 0 ers ' 

description of the input parameters. Uast Sauares - ^ approximating function must be 

1. The maximum generating capacity, G, /, of each utility, is f**r*toc, convex, positive, and increasing. To find the 
i, competing within the market. Note that this capacity is not fnnci } on that best " fits our data whlle satisfying the 
necessarily the total generating capacity of all generators previous constraints, we solve the following quadratic 
owned by a utility. For example, some units may be unavail- program. 

able due to maintenance during certain days. Also, depend- 
ing on the season and availability of fuel, the capacity may 20 ^n, , 
change. The scenario index, s, allows a certain utility to have "^Zj + *if%j + c h ~ 
different possible capacities. The capacity could also change J=1 
depending on t over the planning horizon, usually one week subjct to aj, * o, 
or 168 periods. For example, the capacity over the weekends 

is lower than regular days. 25 ^ * °' 

2. The price functions, f, / at which utility i is willing to xi^a* +xi , *w +c* a 0. 
sell its power at time t under scenario s. The value returned " >J " " 

by f. is assumed to include the operating cost plus the profit 

margin of a utility. Note that f i t s may be given (a) implicitly The previous optimization model can be solved efficiently 

or (b) explicitly. 30 using any quadratic programming approach. See, for 

The Implicit Case example, R. W. Cottle, J. S. Pang, and R. e. Stone, The 

We may not know the exact cost structure of a utility but Linear Complementary Problem, Academic Press, 1992, 

may be able to approximate its operating cost by a quadratic and K. G. Murty, Linear Complementarity, Linear and 

function. For instance, we may observe that a utility, i, Nonlinear Programming, Helderman Verlag, 1988. 

priced its power at time t and under a certain load at X i t s 35 FIG. 2 shows a flow chart of the previous process applied 

$/MWH. By observing the price offered by a utility for its to multiple utilities, 1, . . . ,n, multiple scenarios, 1, . . . ,S, 

power over time, one may construct an approximation of the and multiple periods, 1, . . . ,T. The process starts in function 

cost structure of this utility. The simplest way to do so is to block 201 in which the input data is passed to this procedure, 

collect the different data points of price versus load and fit The input data are the historical prices, X f> 7> , . . . , J,- /, 

a quadratic curve using least squares. Our experience indi- 40 and the corresponding amounts of power generated, x f . 

cates that the quadratic approximation is a good one for most These data are usually available from the power pool and 

utilities. Of course, if one observes a certain profile for the assumed to be stored on the hard disk of our computer. In 

cost structure of a utility, one may consider a different many cases, such data are available through the World Wide 

function to be fit. Web (WWW). At this stage, the counter i is set to one. The 

The problem of "guessing" the price function of a utility 45 process then moves to function block 202 in which the 

is not a simple one in a competitive environment. For counter s is set to one. Function block 203 sets the value of 

example, one may observe a cost of $15 per MWH for a t to one. The process then proceeds to function block 204 in 

certain utility under a 12,000 MWH load on one day while which the quadratic program of equation (1) is constructed 

observing an operating cost of $18 per MWH for the same and solved. We do not discuss the solution approach here, 

utility under a similar load in another day. The reason for this 50 Most optimization packages provide specialized routines for 

discrepancy might be a generating unit that is going through handling such problems. We suggest using the Optimization 

maintenance, a failure of a generator in that utility, a large Subroutine Library of IBM which provides a specialized 

error in the load forecast which resulted in shortage of function, EKKQSLV, for this operation. As a result, we 

generation, or a shortage of certain fuel which may have obtain the values of a^/, b^* and c,.^*, which are used to 

forced the utility to use a more expensive fuel. Regardless of 55 approximate the price function, i.e., production plus profit 

the reason, building the price function f n s of a utility i may margin, for any generation requirement. In function block 

be impossible if one to use a single value for the cost under 205, the value of t is incremented by one. Decision block 

a given load. To handle this difficulty, our model permits the 206 compares t with T. If t^T, the execution is directed to 

use of more than one price function, f^*, at any time period function block 204; otherwise, the execution proceeds to 

t. For a given time period t, we allow the price function to 60 function block 207. In function block 207 the value of the 

be any of the functions f,-/ . . . i t f where S, is the number scenario index s is increased by one. Decision block 208 

of different functions available at time t. That is, we model checks if s is greater than S. If the answer is no, then the 

the cost function of a utility i at each time period t by a set execution is directed back to decision block 203. Otherwise, 

of cost functions. The index s is read as scenario. One hopes the execution continues to function block 209. In function 

that using a set of functions reflects accurately the uncer- 65 block 209, the value of i is incremented by one. Decision 

tainty an observer may have regarding the operating cost of block 210 directs the flow to function block 202 if i is less 

a utility. For example, one scenario may represent a normal than or equal to n; otherwise, the process terminates. As a 
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result, we obtain a set of quadratic functions that can be used 
to approximate the price functions under different scenarios 
and during different time periods. 
The Explicit Case 

The second case we discuss is that in which the price 
function is known explicitly. That is, for each generating unit 
owned by utility i, we know the operating range, the cost 
function within that range, the start-up cost, and any other 
characteristics relevant to a generator. We may also know the 
status of each generator. That is, we know if a generator is 
scheduled for maintenance and if it is unavailable due to 
mechanical problems. Furthermore, one may have historical 
information regarding the probability of failure of each unit. 
All of the previous information implicitly represents the cost 
structure of a utility. That is, we do not have an explicit 
representation of f ir s but we can compute the value of f, v f 
given the load and conditions of the generators. Note that in 
the current electric-power industry, the previous data are 
known to the general public. It is available on the World 
Wide Web and can be obtained in the form of a database 
stored on a floppy disk. However, due to competition, this 
information will not be available in the near future in which 
case we need to rely on the implicit representation men- 
tioned above. 

To build a set of functions, f, that approximate the cost 
structure of a utility under different circumstances, we sort 
the generators in the increasing order of their average 
generating cost. This is done by adding the start-up cost of 
a unit to its operating cost at the maximum operating level 
and dividing by its maximum operating level. To clarify the 
previous approximation, the cost function of FIG. 3 A is 
used. To approximate this function, the start-up cost is added 
to the value of the function at the maximum operating level, 
G. Then, the average cost of this unit is computed by 
dividing the new cost at G by G. As a result, a generator has 
a linear cost function that passes through the origin and its 
operating range is anywhere between 0 and G (see FIG. 3B). 
The generating units are sorted in the increasing order of 
their average costs which results in a convex piece wise 
linear function for operating the system of a utility. FIG. 4 
shows an example of such a function. The first, i.e., 
cheapest, generator is used to cover the load between 0 and 
its maximum capacity, G v The second most economic 
generator is used to meet the electric load between the 
additional load up to the capacity of the second generator. 
That is, it covers the load between G a and G x +G 2 . The 
process is repeated until all the generating units of a utility 
are incorporated into the system. Note that we can build 
different functions to reflect the different operating condi- 
tions of a utility. For example, if some generators are not 
available, then we have a different function. Another 
example is the case in which a generator uses gas as its 
primary fuel. In this case, the cost function can be affected 
by gas prices in the market. Hence, we can create different 
functions (and different approximations) to model different 
market conditions. 

Note that the resulting function is an approximation of the 
cost function of a utility and not its price function. To obtain 
the price function from the cost function we need to estimate 
the profit margin of a utility. The goal of any deregulated 
market is to force the competing entities to reduce their cost. 
The minimum cost at which a utility can sell its production, 
without losing money, is the production cost of its electricity. 
However, this pricing scheme does not generate any profit to 
the utility. To avoid the problem of zero profit, each com- 
petitor usually sets a profit margin that is acceptable for him 
or her. This profit margin may be different for the different 
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periods, t, of the planning horizon and under the various 
scenarios. For example, a profit margin of 1% may be 
acceptable in periods in which the market demand is low 
while a higher margin might be desirable if the demand is 

5 high. The margin can also change with the scenarios s. For 
example, if fuel prices are high, one may not be able to 
enforce high profit margins. 

The profit margin information is assumed to be stored on 
the hard disk of a computer so that it can be accessed by the 

1Q process of this invention when needed. The price function is 
then obtained by adding a linear term that takes into con- 
sideration the profit margin resulting in f t -/. 

To reflect the importance of each scenario and to make the 
model more complete, a probability is associated with each 
scenario s. That is, a probability, p/, is associated with each 

35 price function that reflects the likelihood of its occurrence in 
the future. The value of p/ can be approximated using the 
frequency of each of the functions f t> \ . . . , f i t Sl from the 
past. Of course, the user of the tool according to this 
invention may adjust these probabilities using his or her 

20 intuition if there is a need to do so. These probabilities in 
addition to the parameters representing the functions f ( -/ are 
assumed to be stored on the computer hard disk. 

3. Using the weather forecast, create a load forecast for 
each utility of interest. Usually, these tools require historical 

25 data for the load on the system of a utility and the corre- 
sponding weather conditions, i.e., maximum and minimum 
temperature during the day. Table 2 shows an example of 
this data. 



TABLE 2 

Load vs. Weather Data 



Date 


Max 


Min 


1 


2 


3 


5/25/95 


65 


43 


2101 


2090 


1045 


5/26/95 


68 


44 


2173 


2120 


1932 


5/27/95 


74 


53 


2318 


2292 


2056 



The first column is the date, second column is the maximum 

40 temperature over the twenty-four hours of that day, and the 
third column is the minimum temperature during that day. 
The next twenty-four columns (only three are shown) pro- 
vide the electric load in MWH at each hour of the day. By 
passing this data to a forecasting tool, it can create a model 

45 that correlates the maximum and minimum temperatures to 
the load. These tools are also capable of incorporating the 
effect of time of day and day of year into the forecast. The 
historical data of any utility can be obtained in an electronic 
form from the power pool. In the future, the electric load of 

50 a delivery point, i.e., geographic region, will be posted on 
the World Wide Web. We do not discuss the forecasting 
process in this document. See, for example, R. H. Shumway, 
Applied Statistical Time Series Analysis, Prentice-Hall 
(1988), and A. C. Harvey, Forecasting, Structural Time 

55 Series Models and the Kalman Filter, Cambridge (1991). 
Note that the forecast is given usually for each hour of the 
coming week. The forecasting tool also provides a proba- 
bilistic distribution, i.e., mean and standard deviation, for 
the load at each hour. That is, using the weather forecast, one 

60 can obtain for each utility a probabilistic distribution for its 
load at each hour of the planning horizon. If the forecasting 
tool provides a continuous distribution for the load at each 
time period, the continuous distribution can be approxi- 
mated by a discrete one. The number of discrete realizations 

65 at time t is assumed to be K,. Each demand realization, d, *, 
i-1, . . . ,n, has a probability of it * that reflects the likelihood 
of its occurrence. 
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The weather forecast is usually available on the World total power flowing into a node must be equal to total power 

Wide Web and can be imported by our system easily to pass flowing out of this node. Mathematically, this condition is 

it on to the forecasting system. written as 

4, Transmission cost, c (tV) /, in the different power lines of 

the electrical system. Our notation permits changing the 5 „ „ (3) 

transmission cost from one period to another. It also allows £ y# + = £ yij + 4, i = 1 n 

using multiple cost scenarios depending on the status of the > +1 /=' 
electric system. We assume that for each transmission line 

(ij) we know the maximum transmission capacity MWH whefe d fa lQ be ^ for nodes ^ dQ ^ haye 

which we denote by u ( , vv . If a power line is not available 10 electdc load aQ(J maximum generation, G„ is zero for nodes 

during a certain penod, we can set its capacity, u (lVV , to that do not haye any generation , To satisfy the transmission 

zero. Note that we define a power line by its ends or nodes capacity, we require that 0^y f -^u,. i»l, . . . ,n, j-1, . . . ,n. 

i and j. To simplify the task of representing the electric Also, given the maximum generating capacity of utility i, we 

network and its delivery and generating points, we assume must satisfy me constraint O^x^G,-, i-1, ,n. One can 

that there are n nodes in the network. Each node, i, is 15 enforce other constraints such as the maximum capacity of 

assumed to have a generating utility, i.e., a set of generators, a noc j e 

associated with it and an external demand. Clearly, the upper Due t0 convexity, the minimization of (2) can be solved 

capacity of the generating utility, G,/, at node i can be set easily See> for example, M. S. Bazaraa and C. M Shetty, 

to zero to model a demand center. For a generating node, the Nonlinear Programming: Theory and Algorithms, John 

external load, d, v , can be set to zero. For transmittal or 20 wiley & Sons> Inc ., 1979. if t h e price functions, f ( -, are 

transshipment nodes, we set both generation and demand to quadratic, one can use a speedy quadratic programming 

zero - approach as described by R. W. Cottle et al. and K. G. Murty, 

The transmission cost and capacity are posted on the ibid Note tQat by so \ vmg (2) one obtains the amount of 

World Wide Web as required by the Federal Energy Regu- power t0 be generated and transported from each supplier to 

latory Commission. We assume that our system can access 25 each demand pomU 0 ne also obtains the dual variables or 

the transmission costs and capacities for the different regions me Lagrange multipliers associated with the constraints of 

of interest. equation (3). The dual variable associated with delivery 

Solution Approach p 0 j nt j re p re sents the marginal cost of power at point i. That 

The computer implemented tool of this invention uses the ^ it ^ an estimate of the spot-market price at i under the 

probabilistic distributions of the electric load and of the cost 30 given deman d and price functions. 

function of each utility to build probabilistic distributions for We repeat tDe previous process by sampling a new set of 

the spot-market prices and for the electric trades in the price fu nct i ons aQ( j electric loads then resolving the model, 

system. Due to the use of discrete random variables, i.e., ^ ^mpiing process continues until changes in the distri- 

finite S, and K„ one can consider all combinations of s and bmion of spot prices and the decision variables x iV are within 

k at each hour. The problem of cost minimization while 35 a pre-specifled threshold. Another possibility is to repeat the 

satisfying demand and profit margins is then solved for each pr0 cess for a fixed number, L, of samples, 

one of these combinations to produce a spot-market price ^ flow chart of nG 5 describes the sampling process, 

and to compute the trades. The results can be bucketed to It starts m fu act i ori b i oc k 51 in which the time index, t, is set 

form a distribution of the spot prices and trades. Given that to x ^ process then proceeds to function block 52 in 

the number of combinations at each time period is SpcK,, it 40 which the samp i e counter, 1, is set to one. Note that when- 

may be hard to perform all of the previous calculations in a ever we a new ^ple, i.e., price functions, transmis- 

reasonable time. Instead, one can sample from the electric- sion ^d transmission capacities, the value of 1 is 

load and cost distributions to create a sequence of determin- incremented by one (function block 56). In decision block 

istic models. The sampling process (i.e., simulation) is an 53> we check $ we nave reac h e d the maximum number of 

approximation of the exhaustive one described above. The 45 samp i es , L. If the answer is yes, i.e., 1=L, then the process 

sampling process is repeated until a reasonable approxima- pr0 ceeds to decision block 56 in which we move to a new 

tion is reached of the distributions of trades and spot prices time period j n me planning horizon. If 1 i L, then the process 

or until a maximum number, L, of samples is chosen. moves t0 fun Ct i on block 54 in which a new set of price 

Consider the market equilibrium at time t under a sample factions (one for each utility) and electric loads at the 

of price functions f„ i~l, . . . ,n, transmission cost c, v , 50 different nodes of the network is sampled. This step, func- 

i=l, . . . ,n, . . . ,n, and loads d ff i-1, . . . ,n. We drop the (ion 51ock 54j calls tne sampling procedure of FIG. 6 which 

indices t, s and k to simplify notations. Due to the compe- returns a set of price functions, electric loads, and transmis- 

tition resulting from the bidding process, we expect the spot sion capacities and costs. The data are then used to construct 

price to be the minimum cost at which the demand can be lhe pro bl em of minimizing (2) subject to the constraints of 

met while all suppliers satisfy their profit margins. In other 55 equat ion (3) and other system constraints. We do not discuss 

words, one needs to minimize the cost of supplying power, the so i ut j on process here. Any commercial package that 

supports quadratic programming could be used. As an output 

• V ■ V V (2) for the optimization, we obtain the amount of power flowing 

fi{Xi) + 2j 2j c u y v- in the different lines of the network, / ( , ^ and the marginal 

J 60 costs, n?f „ corresponding to equation (3) at the different 

nodes of the system. The process is then directed to function 

The decision variable x f represents the amount of power block 53 in which a new sample is generated. After reaching 

produced by supplier i while y lV represents the amount of the maximum number of samples in decision block 53, the 

power flowing from i to j. Note that if two nodes are not process moves to function block 56 in which the time 

connected directly with a transmission line, then the capac- 65 counter, t, is incremented by 1. In decision block 57, if t^T, 

ity of this line is set to zero. The minimization of (2) must then the process goes back to function block 52; otherwise, 

be solved subject to Kirchhoff-law which requires that the the process terminates. At the end of the process of FIG. 5, 
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we have for each time period a set of L values for the 
spot-market prices at nodes i-1, . . . ,n. We also have a set 
of possible power flows between the different nodes in the 
transmission system. 

Note that function block 54 of FIG. 5 calls the process of 
FIG. 6 which performs the sampling. The sampling process 
is performed using a uniform random number generator. 
Uniform random number generators are available on most 
computers and under any programming language. For 
example, the C programming language provides the function 
rand as part of its standard library header. See B. W. 
Kernighan and D. M. Ritchie, The C Programming 
Language, Prentice Hall, 1988. The probability of the dif- 
ferent scenarios, p/, is used in conjunction with the electric- 
load probabilities, ji,*, to come up with a sample. The 
process of FIG. 6 receives from the calling process (function 
block 54 of FIG. 5) the value of t. In function block 601, a 
random number x between 0 and 1 is sampled from the 
uniform distribution. Function block 602 defines a variable, 
y, and sets its value to zero. The variable y is used to 
compute the cumulative joint distribution at any given k and 
s. The value of the scenario counter, s, is set to 1. Function 
block 603 sets the value of k to 1 . In decision block 604, the 
value of the cumulative probability, y, is compared to x. If 
x is less than or equal to y, then the functions f,- / and loads 
d,-/ are returned to the calling routine (function block 54) to 
be used in the optimization problem. Otherwise, the value of 
k is incremented by 1 in function block 605. Decision block 
606 checks if k has reached its maximum limit of K r If 
k^K„ then the value of the cumulative probability, y, is 
increased by p/n;* -1 and the execution is sent back to 
decision block 604. Otherwise, the value of s is increased by 
1 in function block 608. Function block 609 updates the 
value of y by adding p/ _1 n: * _1 to it. The execution then 
loops back to function block 603. Function block 610 is 
reached from decision block 604 when the value of x is less 
than the cumulative distribution at the point (k,s). 
System's Output 

At termination of the simulation runs, we have the spot 
prices or marginal costs, X 1 1 r , at each point i, i-1, . . . ,n, and 
at each time period of the planning horizon. Here, the 
superscript 1 indicates the sample index. We also obtain the 
power flow, y'^,, in the different transmission lines of the 
power network. Note that if y^,, is greater than zero, then 
node i is transmitting power to node j under sample 1 of our 
simulation. Of course, if y^^O, then /(mv^O* ^° & et 
probabilistic distribution of the spot-market prices, \* i/t we 
perform the following steps. 

1. Create ranges for \ l i r For example, we may want to 
have intervals of S5/MWH width. That is, we have the 
intervals (O-5],(5-10],(lO-15] 

2. Count the number of r that falls within each range. 
We normalize these numbers to obtain a probability for each 
range. 

3. The resulting histogram is an approximation of the 
probability distribution of the marginal cost at delivery point 
i at time t. 

The previous process is called bucketing and is described 
in the flow chart of FIG. 7. The process requires a number 
of buckets (i.e., ranges) which is defined by the user and the 
lower limit, u ; ., of each range. We can assume that there is an 
extra range (after the last one defined by the user) that covers 
the range up to infinity. Function block 71 sets the number 
of observations (so far) in each range or bucket to 0. The 
process then moves to function block 72 in which the value 
of yJ i { is compared against the values of u, to find the index 
j for which Mj^Ki t t= u j*i' Th e value of the counter m ; - is 
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incremented by 1 in function block 73 to account for sample 
j. Function block 74 increases the value of 1 by 1; i.e., we arc 
ready to process a new sample if any. Decision block 75 
checks if 1 is less than L. If this is the case, then we have 

5 more samples to study and the process moves back to 
function block 72. Otherwise, we completed classifying \ l it 
into the appropriate buckets. Finally, function block divides 
the frequency, m ; , by L to obtain the relative frequency for 
each range or the probability. 

We can follow similar steps to create probabilistic distri- 
butions for the amounts of power traded, y 1 ^^ across the 
different power lines, (i,j). The resulting distributions can be 
helpful to the decision maker in deciding the amount of 

is power to be generated and its price. For a power broker, 
these distributions can be used to sign buying and selling 
contracts with different suppliers and consumers. 

While the invention has been described in terms of a 
single preferred embodiment, those skilled in the art will 

20 recognize that the invention can be practiced with modifi- 
cation within the spirit and scope of the appended claims. 

Having thus described our invention, what we claim as 
new and desire to secure by Letters Patent is as follows: 

25 1. A computer power trading and forecasting process for 
forecasting spot-market prices of electric power and trading 
transactions at different delivery points comprising the steps 
of: 

inputting maximum generating capacities of each utility 

30 competing within the market, price functions at which 
a utility is willing to sell its power at a given time 
period, a weather forecast, and a transmission cost and 
transmission capacity in different power lines of the 

35 electrical system; 

using probabilistic distributions of electric load and the 
price function of each utility and transmission cost and 
transmission capacity between different utilities to 
build probabilistic distributions for the spot-market 

40 prices and for electric trades in the system; and 

outputting probabilistic distributions for spot prices and 
trades at a plurality of time periods of a planning 
horizon and for a plurality of delivery points to enable 
a power broker to effectively manage risk in buying and 

45 selling electric power with different suppliers and con- 
sumers. 

2. The computer power trading and forecasting process 
recited in claim 1 wherein the maximum generating capacity 
of each utility changes under different scenarios. 
50 3. The computer power trading and forecasting process 
recited in claim 1 further comprising the step of using the 
input weather forecast to create a probabilistic distribution 
of load forecasts for each utility of interest. 

4. The computer power trading and forecasting process 
55 recited in claim 3 wherein the input price functions are 

derived implicitly from production plus profit margin. 

5. The computer power trading and forecasting process 
recited in claim 4 wherein a price function of a utility is 
approximated by a quadratic function which is convex, 

60 positive and increasing. 

6. The computer power trading and forecasting process 
recited in claim 3 wherein the input price functions are given 
explicitly as a cost function of the generating units. 

7. The computer power trading and forecasting process 
65 recited in claim 3 wherein probabilistic distribution of the 

transmission costs and capacities of power lines are created 
from historical data. 
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8. The computer power trading and forecasting process 
recited in claim 3 wherein the probabilistic distribution for 
the spot-market prices and for electric trades in the system 
are created by sampling from electric-load, price functions 
and cost distributions to create a sequence of deterministic 5 
models. 

9. The computer power trading and forecasting process 
recited in claim 8 further comprising the steps of: 

minimizing a cost of operating an electric utility system 
under each sampled set of price functions, transmission 10 
costs and electric loads; and 

storing a resulting output. 

10. The computer power trading and forecasting process 
recited in claim 9 wherein sampling from electric-load, 
transmission costs and cost distributions is performed using 15 
a random number generator. 
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11. Hie computer power trading and forecasting process 
recited in claim 3 wherein the probabilistic distributions for 
spot prices and trades is performed by a bucketing process. 

12, The computer power trading and forecasting process 
recited in claim 11 wherein the bucketing process comprises 
the steps of: 

creating ranges of spot prices; 

counting a number of spot prices that fall within each 
range and normalizing the counted numbers to obtain a 
probability for each range; and 

generating a histogram of the probabilities of each range 
as an approximation of a probability distribution of the 
marginal cost at each delivery point at a given time. 

***** 
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